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The WA98 collaboration in the CERN SPS has reported an excess of photons over those 
originating from the decays of final hadrons in the lead-lead collisions. These photons 
can originate either from primary interactions of partons from colliding nuclei or from 
secondary interactions among produced particles. Photons produced in the secondary 
interactions, often called thermal photons, can be calculated by using thermal production 
rates and equilibrium hydrodynamics for the evolution of the expanding matter. I will 
review the main features of hydrodynamical studies for the WA98 data. The data can 
be reproduced both with or without a phase transition to the QGP, but high initial 
temperature, over the values predicted for the phase transition temperature, is required 
by the data. I will also show a prediction for the photon excess for central gold-gold 
collisions at the Brookhaven RHIC collider. In this prediction, the initial state for the 
hydrodynamical expansion is obtained from a perturbative QCD calculation. 

1. Introduction and Theoretical Framework 

Photons are suggested to be an ideal probe for the hot initial state of the relativistic 
heavy ion collisions, since, due to the large mean free path in a hadronic matter, they 
escape from the collision zone right after being emitted without any further interactions. 
Unfortunately, photons emitted at early stages of the collision are difficult to obtain ex- 
perimentally, since the decays of the final state hadrons, especially neutral pions and etas, 
dominate the production of photons in the collision. So far only the WA98 collaboration 
at the CERN SPS has reported an excess of photons over those originating from decays of 
the final state hadrons [[[],§. Photons are measured also by the PHENIX collaboration 
at Brookhaven RHIC collider, but more statistics is needed to reach conclusions [ §. 

In a hydrodynamical description of relativistic heavy ion collision the secondary inter- 
actions are assumed to lead to thermalization of the fireball. With this assumption, the 
evolution of the system can be described by equilibrium hydrodynamics and transverse 
momentum spectrum of the thermal photons emitted during the expansion can be calcu- 
lated if thermal emission rate in the matter is known. Experimentally measured direct 
photon spectrum has a contribution from primary interactions as well [ |J, referred as 
pQCD (or prompt) photons in following. 

Several groups [ §, 0, @, H have performed hydrodynamical studies to explain WA98 
data. The authors of [ [5] use a (2+l)-dimensional hydrodynamical code that takes into 
account also the finite longitudinal size of the initial collision zone, whereas the studies in [ 
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H) H? El are based on the boost-invariant approximation. In the boost-invariant scenario 
very hot initial state [ [| , considerable initial radial velocity [ [7], ^], [J or modification of 
hadron masses [ [| is needed to explain the WA98 data. I have made a choice to omit 
discussion on modifications of hadron masses [ |9[] and also chosen references [ || [7j to 
represent boost-invariant studies in this note. 

Important improvements have been obtained in a perturbative calculation of the ther- 
mal emission rate of photons from the thermalized QGP [ [H| [11], In hadron gas 
mesonic processes tttt — > p7 and np — > pj, described with a pseudo-vector Lagrangian, 
are included in the photon emission rate [ [13|] . The up scattering channel gets a large con- 
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tions of the emission rates are provided by the authors of [ [TJ], [15 
are not included in the rates. 



PI 

The effects of baryons 



2. Details of Hydrodynamics and Results 

In section ^TT| results of three different groups [ are compared to the WA98 data 
[ H' H ■ Groups [ |6|, |7j use 2-loop photon emission rates in QGP [ [16], [L7j while in [ U 
resummed rate [ [121 , complete in order a s , is used. In works [ |], || contribution of the 
pQCD photons is taken from [ [Ti|]. The authors of [ [7J present their own estimate for the 
pQCD photons. 

In section a prediction for direct photon spectrum in Au+Au collisions with y/s = 
200 AGeV is given. Results for hadronic observables in this study, where the initial state 
is based on perturbative QCD + saturation model, can be found in [ [T^, |21J. For this 
energy, hard pQCD photons from the primary interactions are studied in [ 

The following discussion will concentrate on the role of the initial state in hydrodynam- 
ical models, since especially high transverse momentum k t part of the thermal photon 
spectrum is sensitive to the initial conditions. In particular, temperature T in the initial 
state plays a big role due to the Boltzmann suppression e~ kt l T in the emission rates. 



2.1. Comparison with the WA98 Data 

In boost-invariant hydrodynamical studies one must give the initial energy density 
e(r, To) and velocity v (r, To) distributions in the transverse plane at fixed formation (ther- 
malization) time To, that are the initial conditions for the hydrodynamical evolution. Let 
us start with a case v(r, tq) = for all r. For the shape of the initial energy density 
e(r, To) various choices can be found in the literature. For example, it can be proportional 
to wounded nucleon distribution [ |6|, Woods-Saxon distribution [ [?J or nuclear overlap 
function [ [0J. Within these choices, the formation time r remains as the main uncer- 
tainty in the boost-invariant picture. At SPS energies, r is expected to be ~ 1 fm/c 
based on the geometrical argument, that it takes a time ~ 2Ra/j, where Ra is the nu- 
clear radius and 7 the Lorentz gamma-factor, for the colliding nuclei to pass through each 
other. Fixing To = 1 fm/c leaves only the normalization of the initial energy density as a 
free parameter, that can be fixed from experimentally measured hadron spectra. 

Studies by the authors of [ |(| reveal a problem in this simple approach. If initial time 
is chosen to be of the order 1 fm/c, the initial state is way too cool to reproduce high 
transverse momentum, k t , part of the spectrum as measured by the WA98 collaboration. 
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Figure 1. Direct photon production in ref. 
[ § with the WA98 data. QM refers to 
emission in QGP and QGP part of mixed 
phase, HM likewise for hadronic matter. 
T c is a critical temperature in EoS, r 
thermalization time and T average initial 
temperature. 
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Figure 2. Direct photon yield calculated 
for EoS with phase transition and linear 
initial velocity distribution in ref. [ 0] with 
the WA98 data. Initial parameters for dif- 
ferent curves in the figure are presented in 
the table |l|. Estimate for pQCD (prompt) 
photons is included to the curves. 



Since decreasing of initial time considerably raises the temperature in the initial state| and 
final state hadron spectra are fairly insensitive to small changes of To, the data seems to 
suggest very rapid thermalization and hot initial state within the boost-invariant scenario. 

Figure [l] shows the main results in [ . The data is well reproduced with initial time 
r = 0.2 fm/c which corresponds to an average initial temperature (T) ini = 335 MeV in 
their choice of the initial state. This temperature is considerably higher than the expected 



critical temperature T c = 150 — 180 MeV from recent lattice calculations [ [22] and hence 
the authors of [ || conclude that this confirms the formation of the QGP. However, in my 
opinion, one should consider the initial time tq = 0.2 fm/c carefully, although the WA98 
data seem to favor it. The geometrical argument above, on the transit time of the nuclei, 
suggests that at such a short initial time primary collisions are still taking place, and only 
part of the matter has been produced [ ^] . 

Let us turn to discuss effects from non-zero initial radial velocity. The authors of [ 
[7| have studied the effects of radially increasing initial velocity t>(r, r ) = (r / Ra)Q(R — 
r)y max , where V max is treated as a free parameter. They have performed a simultaneous 
X 2 analysis of 7r° and 7 spectra measured by the WA98 collaboration [ |l| to find the best 
combination of the thermalization time Tq and the velocity parameter V ma/X . Their results 
[ [7J with a phase transition to QGP are presented in figure |2|. Curves in the figure |2| 
correspond to the choices of the initial parameters presented in the table [I]. The results 
show clear interplay between the initial temperature and the strength of the initial radial 
velocity: thermal photons with large kt are emitted at first moments of the evolution when 
the temperature is high. Lowering of initial temperature, or equivalently increasing the 
thermalization time Tq, requires an introduction of initial radial velocity, which boosts the 



2 If radial expansion is omitted, T(r, r) = T(r, r ) x (t/tq) 1 ^ 3 for massless ideal gas. 
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Table 1 

Choice of the initial parameters in [ [7]]. 



(T) ini [MeV] 


r [fm/c] 






Curve in the figure ^| 


180 


2.0 


0.40 


solid 


200 


1.5 


0.38 


dashed 


230 


1.0 


0.29 


dotted 


250 


0.7 


0.26 


dash-dotted 



emitted photons to larger k t . With this compensation, initial temperature can be lowered 
below 200 MeV and the thermalization time rises to values 7o>l fm/c. These values of 
the temperature are quite close to the critical temperature and hence the authors of [ [7j 
conclude, that, considering uncertainties of the model, formation of QGP or its absence 
cannot be distinguished. 

The origin and the magnitude of the initial radial velocity is a somewhat uncertain. The 
authors of [ [F| explain its origin "as a result of significant radial energy gradient coming 
from the shapes of the colliding nuclei" . They do not present a (quantitative) model that 
could explain how this leads to a collective motion into the radial direction at initial time 
t . In this sense, I would consider y max to be an extra fitting parameter in the model. 
In any case, initial radial velocity can not be large. The authors of [ |5j have made a 
boost-invariant calculation with initial radial velocity profile taken from [[?]. They found 
a good agreement with the photon data by choosing initial parameters To = 1.0 fm/c 
and ^/ max = 0.30, but the resulting tt° spectrum clearly overshoots the experimental one. 
It should be noted that the initial energy density profile — normalization is fixed from 
multiplicity in both cases — and the equation of state are different in these two studies, 
and in [ [j| chemical and kinetic freeze-out temperatures are the same. These differences 
may explain the discrepancy on the 7r° spectrum. 

The authors of [|5j used a (2+l)-dimensional hydrodynamical code, where also the lon- 
gitudinal extend of the fireball is finite and the longitudinal velocity evolves dynamically 



from a given initial state. Details of this model can be found in references [ ^4|, |25|, pq| . 
My first remark is, that there is no uniquely defined initial time, when the assumption 
of boost-invariant flow is relaxed. In the boost-invariant case the source is pointlike in 
the zt-plane and one can trace all the particles back to the collision point. This cannot 
be done when the longitudinal flow profile of the system is not that of the scaling flow, 
v z = z/t. To estimate the thermalization time one can relate a timescale to the longitu- 
dinal extent by considering the collision geometry. In the model [ ^5j the initial length of 
the system was chosen to be 3.4 fm. Since 7 ~ 10 at SPS energies, the system is roughly 
1.3 fm thick when the colliding nuclei overlap and it would take ~ 1 fm/c for the system 
to reach longitudinal extent of 3.4 fm [ [27j]. The rapidity distribution of initial energy den- 
sity is assumed to be Gaussian. The energy per unit transverse area, e(r) = / dzT 00 (r, z), 
is assumed to equal the energy per unit transverse area of the incoming nucleons. To 
convert this distribution to spatial energy density, one must define initial velocity profiles. 
Initial conditions, chosen to reproduce the observed hadron spectra, depend also on the 
equation of state (EoS). In [ [5] EoS H describes a hadron resonance gas without phase 
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Figure 3. Number density dN/dy of neg- 
atively charged hadrons function of 
rapidity y with the NA49 data [ |, ||]. 
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Figure 4. Thermal photon emission with 
the contribution of pQCD photons (upper 
set, scaled by factor of 100) and without 
(lower set) in ref. [ |J. Solid line show the 
results with QGP formation and dashed 
line in purely hadronic scenario. 



transition and EoS A has a first order phase transition to QGP. Initial radial velocity is 
set to zero. For the longitudinal velocity profile two different choices were studied [ ^BJ; 
in the case called IS 1 the longitudinal flow rapidity is assumed to increase linearly with 
the distance z. This gives a large peak, by a factor of two over the Bjorken estimate, to 
the initial energy density. In the other choice, IS 2, the z dependence is nonlinear and 
the resulting energy density profile is almost flat. 

Combinations IS 1 + EoS A and IS 2 + EoS H reproduce simultaneously the NA49 data 



for the rapidity distributions of hadrons [ ^8f, and the WA98 data for neutral pions and 
direct photons [ [L| . Figure |3| shows the results for measured rapidity density of negative 
charged hadrons, comparison to 7r° spectrum can be found in [ |5| and figure |] shows the 
results for direct photons. Maximum and average temperatures in z = in the initial 
state are given in the table |2|. 



Table 2 

Maximum and average initial temperature in [ [5[]. 





IS 1 + EoS A 


IS 2 + EoS H 


T max [MeV] 
(T) ini in z = [MeV] 


325 
255 


245 
213 



At the given temperature, the thermal emission rate of photons is larger in hadronic 
matter, which makes IS 2 with lower initial temperature more suitable for purely hadronic 
scenario. Choosing IS 1 + EoS A leads to a high initial temperature, but gives a good 
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agreement also with the high kt part of the measured photon spectrum, as shown in the 
figure [|. Despite the ~ 50 MeV difference in the initial temperature, later evolution is 
so similar in these two choices for the initial state, that hadronic spectra are equally well 
reproduced. These observations reflect the ambiguity in the choice of the initial state. 



After pQCD photon spectrum from the ref. [ |18| is added to the thermal one (upper 
set of curves, scaled by a factor of 100, in the figure f|), the authors of [ |5| find that the 
WA98 data can be reproduced both with or without a phase transition to QGP. However, 
the maximum (average) temperature 245 (213) MeV in the initial state with hadron gas 
is fairly high, and clearly over the expected phase transition temperature 150—180 MeV 
[ ^2f . Also the shape of the solid curve in the figure £|, which corresponds to the results 
from IS 1 + EoS A with the phase transition, seems to follow the data better. 

2.2. Thermal Photons at RHIC Collider 

References [ [111 EU| show results for hadronic observables from boost-invariant hydro- 



dynamics, when the initial state is obtained from a pQCD + saturation model. With this 
model for the initial state, only free parameter in the hydrodynamical description is the 
decoupling temperature. The initial time and energy density, including normalization, 
can be calculated from the model when the center of mass energy of the collision, ^/s, 
mass number of the colliding nuclei, A, and the experimental centrality cut are given. 
Construction of the initial state is explained in [ [OJ, where also results for global ob- 
servables (e.g. multiplicity) are shown. In [ ^U] the effect of the decoupling temperature 
is discussed and results are compared to hadron spectra measured by the PHENIX col- 
laboration [ pfflfl . Following results for thermal photons are calculated with initial 
parameters yfs = 200 AGeV and A = 197, which lead to r = 0.17 fm/c, a(E T ) = 86.33 
mbGeV and A efi = 178 for a 6 % most central collisions. For details see [ |H| . 

The initial time is determined from the saturation momentum scale, p sat , of the pertur- 
bative minijet calculation [ |31~1| . This leads to a very short initial time r = l/p S at = 0.17 
fm/c compared to values >0.6 fm/c, that are used in many recent hydrodynamical studies 
at RHIC energies, for a review see e.g. [ |32| . The short initial time r = 0.17 fm/c is 



consistent with the geometrical argument presented in section 27T . The produced minijet 
system also looks thermal from the point of view of the number of the partons and the 
energy per particle at the initial time [ [HJ. Starting hydrodynamics at tq = 0.17 fm/c can 
be considered to give an upper limit for the initial temperature and the thermal photon 
production. 

The solid line in the figure [5] shows the transverse momentum spectrum of thermal 
photons. The dashed line shows contribution from quark-gluon plasma and QGP part of 
the mixed phase, and the dotted line likewise from hadron gas. The plasma contribution 
dominates the thermal emission for k t > 3 GeV. Comparing results from SPS and RHIC 
energies one can see roughly a factor of 100 increase in thermal emission at kt = 4 GeV, 
but the difference is less than order of magnitude at kt = 1 GeV, which reflects a clear 
increase of the maximum temperature. 

To analyze the effect from initial temperature (thermalization time) I have scaled the 
initial state from tq = 0.17 to 1.0 fm/c with longitudinal Bjorken expansion so that 
entropy in unit rapidity does not change. Neglecting transverse expansion before r = 1.0 
fm/c is not physical, but in this way one can define a cooler initial state in such a way 
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Figure 5. Thermal photon spectrum in 
Au+Au collision with y/s = 200 AGeV. 
Plasma refers to the contribution from 
QGP and QGP part of the mixed phase. 
Also hadron gas includes emission from 
hadronic part of the mixed phase. 



Figure 6. Effect of initial time tq to 
the thermal photon emission. Spectra of 
pQCD photons, with two different values 
for intrinsic transverse momentum (pf), 
are from work by A. Dumitru et al. j 
[|. Initial temperatures are in the ta- 
ble |. 



that the results for hadronic observables do not change much. Table |3] shows the changes 
in the maximum and average values of the initial temperature and energy density. Effect 
on the thermal photon spectrum is presented in the figure |6|. 

Increasing initial time from r = 0.17 to 1.0 fm/c reduces the thermal emission by a 
factor of 20 at k t = 4 GeV. Unfortunately, the amount of thermal photons in this region 
may be hard to resolve experimentally, because pQCD photons from primary interactions 
may dominate the direct photon spectrum there. The analysis of the pQCD photons can 
be found in [ ^TJ. In [ [2l| the centrality cut is 10 % instead of 6, so the pQCD results 



presented in the figure ^| are multiplied with a factor of 1.1 [ [£J. Results for two different 
values of the intrinsic transverse momentum, (p 2 ) = and 2.4 GeV 2 , are given. 

These theoretical results on excess photons at RHIC can be summarized as an upper 
limit given by the sum of the solid and dashed lines and a lower limit given by the sum 
of the dotted and dash-dotted lines. Uncertainties are fairly large, but in any case one 
should see a change in the slope in the region k t ~ 2.0 — 3.5 GeV, above which the pQCD 
photons start to dominate the spectrum. In principle a very clean data around k t ~ 2 
GeV could fix the value of the initial time, but one should keep in mind the possible 
effects from initial transverse velocity and longitudinal geometry discussed in the section 
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Table 3 



Maximum and average initial temperatures and energy densities in Au+Au collisions with 
y/s = 200 AGeV. 





T max [MeV] 


(T) ini [MeV] 


e m ax [GeV/fm 3 ] 


(e) ini [GeV/fm 3 ] 


r = 0.17 fm/c 


580 


340 


208 


60 


t = 1.0 fm/c 


320 


220 


20 


6.8 



27T| . These effects may be smaller at RHIC energies, because rapid thermalization means 
shorter time interval to build up the initial transverse velocity. Also the rapidity range is 
broader at RHIC, and hence scaling hydrodynamics should work better. 



3. Conclusions 

I have reviewed hydrodynamical studies to explain the WA98 excess photon spectrum. 
The role of the initial state and the uncertainty in the values of initial temperature were 
emphasized in the discussion. The data can be explained with or without QGP phase 
transition, but high initial temperatures, well above the predicted critical temperature, 
are favored. In the case without boost-invariance the shape of the spectrum is closer to 
the data, when a phase transition is assumed. Still, both theoretical and experimental 
uncertainties are somewhat too large to rule out either alternative. 

Using boost-invariant hydrodynamics with an initial state calculated from a pQCD + 
saturation model, theoretical estimates of the upper and lower limits for the thermal pho- 
tons were given and compared with the results for the pQCD photons [ [|, ^TJ. 
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